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Abstract 



One of the objectives of the vast physics programme of the recently comissioned WASA- 
at-COSY facility is the study of fundamental symmetries via the measurements of the 77 
and 7/ mesons decays. Especially interesting are isospin violating hadronic precesses of 
these mesons into 3tt systems driven by the term of QCD Lagrangian which depends on 
the mass difference of the u and d quarks. 

When an r? or an 7/ meson is created in the hadronic reaction signals from such de- 
cays may be significantly obscured by the prompt production of tt mesons. In this thesis 
we present the estimation of the upper limit of the background due to prompt pion pro- 
duction for the 7/ — > 3ir° and 7/ — > ir + ir~ir° decays. Using the data from proton-proton 
collisions measured by the COSY- 11 group we have extracted differential cross sections for 
the multimeson production with the invariant mass corresponding to the mass of the 7/ 
meson. Based on these results and on parametrizations of the total cross sections for the 
7/ meson as well as parametrization of the upper limit for the prompt tt + tt~tt° production 
in the collisions of protons we discuss in details the feasibility of a measurement of the r/ 
meson decay into 3ir channels with the WASA-at-COSY facility. Based on the chiral uni- 
tary approach the value of the branching ratio BR(r/ — > 7r + 7r~7r°) was recently predicted 
to be about 1%. We show that the WASA-at-COSY has a potential to verify this result 
empirically. 

Furthermore, we discuss the possible usage of the time signals and of the energy loss 
measurement in the forward part of the WASA-at-COSY detector for the determination of 
the energy of the forward scattered charged particles. We briefly describe the detectors to 
be used for this purpose and introduce the computational algorithm which was developed 
to be applied with this technique. 



"There is no expedient 

to which a man will not go 

to avoid the labor of thinking." 



"Nie ma takiego fortelu, 

do ktorego nie odwolaiby sie^ czlowiek, 

aby uniknac pracy zwanej mysleniem." 



Thomas Alva Edison (1847 - 1931) 
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1. Introduction 



The Standard Model is nowadays a well known and established theory describing the 
particles and theirs strong, electromagnetic and weak interactions. Therefore, tests of the 
applicability of this model are very important and are carried out in many particle physics 
laboratories. In particular particle physicists endeavour to discover phenomena which the 
Standard Model does not describe. One of such laboratory is WASA-at-COSY [1] where 
we investigate the limits of applicability of the Standard Model by studying production 
and decays of pseudoscalar mesons like 77 and 77'. The examination of these meson produc- 
tion and their decays give us a chance to probe fundamental symmetries like C (charge 
conjugation), P (parity), T (time reversal) and their combinations and to determine the 
parameters of the Standard Model like for example the quark masses. 

The masses of the light quarks (up -u, down - d, strange - s) are known only approx- 
imately due to the fact that quarks are confined inside hadrons and cannot be observed 
directly as physical particles. But there are several indirect possibilities of establishing 
these masses by measuring the mass differences and mass ratios. The mass difference of 
light quarks may be derived from the partial widths of the decay of the 77 and 77' mesons into 
7T + 7r~7r and theirs total widths. The possibility of determining of the difference between 
the mass of the light quarks motivates us for studing these decays. 

Detailed studies of 77 — ► tt + it~it decay was conducted by the WASA/CELSIUS collab- 
oration [2], but so far the decay 77' — ► 7r + 7r~7r° has never been observed. In both of these 
processes the isospin conservation is violated [3]. But despite this fact in case of the 77 
meson the branching ratio for the 77 — ► tttttt decay is in the order of 50%. However, for the 
77' meson the situation is quite different the branching ratio for the 7r 7r°7r° is at the permil 
level and for the 7r + 7r _ 7r° system so far only an upper limit of 5% was established [4]. 

For the decay of both 77 and 77' mesons there exists a physical background from the 
direct production of three pions via the pp — ► ppir + iT~iT reaction channel. In case of the 
77 meson the signal to background ratio is large, amounting to about 10 for tagging by 
means of the missing mass technique with a resolution of a few MeV, and permits a clear 
identification of the 77 — > 7r + 7r~7r° decay. However, this ratio is expected to be worse by 
more than three orders of magnitude for the 77' meson making the investigations much 
more challenging experimentally, especially since the hadronic production cross section is 
by about a factor of thirty smaller for the 77' meson in comparison to the 77 (a v « 30ay) 
meson at the same excess energy. 

The main aim of this thesis is the determination of the optimum beam momentum 
for the measurement of the branching ratio for the 77' — > 7r + 7r _ 7r° decay channel via the 
pp —> pprj' — > ppir + ir~7r — ► ppir + Tr~ / yj reaction chain. The measurement is planned to be 
carried out in the near future and its result will be used for the derivation of the light quarks 
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mass difference (m^ — m^). The experiment will take place at the Research Centre Jiilich in 
Germany and will be conducted by the WASA-at-COSY collaboration. The identification 
of the reaction is based on selecting events where the tt + tt~tt° system was produced. For 
this subsample of data the fractions corresponding to the direct 7r + ir~ir° production and to 
the decay of the rf meson will be extracted based on the missing mass of the two outgoing 
protons registered in the Forward Detector of the WASA-at-COSY system. 

To calculate the number of expected rf events we have parametrized the total cross 
section for the rf meson production and have etablished a parametrization for an upper 
limit of the background production. Further on we have parametrized a missing mass 
resolution taking into account effects which are related with the energy resolution of the 
forward detector and the beam and target spread. In the case of the rf meson the missing 
mass resolution is at present not satisfactory for studies of rare decays due to the high 
energies of the protons. Therefore, we have proposed the usage of the time-of-flight method 
(TOF) for the identification of charged particles emitted in forward direction based on time 
signals from the scintillator detectors. 

In the next chapter of this work we will outline difficulties of the mass determination for 
light quarks. We will quote the idea of an indirect determination of quark masses proposed 
by D. Gross , B. Treiman and F. Wilczek [5], and by H. Leutwyler [6]. 

The measurment and identification method of the decay rf — > 7r + 7r~7r° as well as the 
WASA-at-COSY detector setup will be presented in Chapter 3. 

The first section of Chapter 4 describes how the total cross section for the production 
of the rf meson depends on the excess energy near the kinematical threshold. Further on in 
this chapter the method of estimating an upper limit of the background under the rf peak 
in the missing mass distribution will be given [7]. Next, applying the parametrization of 
the energy resolution of the detector, the accuracy of the branching ratio determination as 
a function of excess energy and time will be shown. For the calculation we consider a range 
of values for the upper limit of 5% down to a value lower by one order of magniude (0.5%). 

In chapter 5 we will propose a method of particle energy reconstruction by measuring 
the time in the Forward Detector. For the purpose of this thesis we will consider only thin 
scintillators, but in the future the reconstruction will be based on time signals form all 
detectors. First we will show results of simulations for an energy resolution based on TOF 
measurements and describe how passive material of the detector influences the precision 
of the TOF determination. As a result of the simulations the missing mass distributions 
calculated based on the TOF measurment will be shown. In the last section of this chapter 
we will describe the computational algorithm which was developed to reconstruct the 
energy of particles based on the time signals from scintillators and energy losses in the five 
independent layers of the range hodoscope. 

Chapter 6 summarises the whole thesis and brings the conclusions and remarks regard- 
ing the beam momentum and the duration of the planned measurement of the rf — > 7r + 7r~7r° 
decay. 
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The work is supplemented with Appendics where section [A] introduces basic informa- 
tions about the rf meson and the SU(3) symmetry. Appendix [B] describes the missing mass 
technique used for tagging the rj meson production. Section O provides an analytical calcu- 
lation of the fractional momentum resolution as a function of the fractional time resolution. 
In appendix iDl the parametrization of the total cross section for the 7r + 7r~7r° production 
will be shown. Section [E] introduces the parametrization of the Final State Interaction 
(FSI) and production dynamic of the multipion system. The next section presents simula- 
tions of the 37r production in proton-proton collisions using models described in section lEl 
The last appendix shows schemes of programs used for the simulations. 
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2. Relation between the partial width T^^+^^o 
and the u and d quark mass 

The determination of the light quark masses is one of the important goals of hadron 
physics experiments, and we intend to contribute to their estimations by determining the 
quark mass difference — m u which induces an isospin breaking. Studing the isospin- 
violating decays 77(7/) — ► 7r + 7r~7r° and 77(7/) — > 7r 7r°7r° was pointed out as an accurate way 
of extracting the quark mass difference [5,6]. 

The quark masses are one of the Standard Model parameters and their values depend 
on how they are defined. The field theory which describes the strong interaction between 
gluons and quarks is the Quantum Chromodynamics (QCD). A general QCD Lagrangian 
for N flavors reads: 

N 1 
Lqcd = Yl ^{iD - m k )q k - -G^G^, (2.1) 

fc=i 

where ra k denotes the quark masses, D indicates the gauge convariant derivative and G^ 
represents the gluon field strength. 

In the low energy regime where the strong coupling constant a s is large, the perturbative 
approach of QCD cannot be applied in the same way as for high energies. Thus other 
methods must be applied in low energy hadron physics. One of these approaches is an 
effective field theory of QCD at low energies - the Chiral Perturbation Theory (ChPT). 
It is based on the observation that, in the low energy region, the relevant and effective 
degress of freedom of strong interactions are hadrons composed of confined quarks and 
gluons. This leads to the effective Lagrangian which is formulated in terms of the effective 
degrees of freedom: 

^QCD = + (2.2) 

We can use this effective Lagrangian instead of the formula (|2.1|) . The Lo term is the part 
of flavor symmetry of QCD and L m contains a contribution: 

- \( m d ~ m u )(uu - dd) (2.3) 

responsible for the isospin changing in QCD (AI = 1). One can identify the degrees of 
freedom as eight Goldstone bosons which are n, K and r\ mesons from the pseudoscalar 
meson nonet. This effective Lagrangian breakes down spontaneously the chiral symmetry 
SU(3)l x SU(3)r to SU(3)v- It is believed that this violation of the chiral symmetry 
constitutes the source of the mass generation in QCD. The standard ChPT provides an 
accurate description of the strong and electroweak interactions of the pseudoscalar mesons 
at low energies. 
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and the u and d quark mass 



The rj and rj mesons are members of the pseudoscalar nonet, and play an important 
role in the understanding of the low energy QCD. Due to the mixing of the rj and rj' fields 
the 77 treatment in the standard ChPT is complicated. The rj' meson is related with the 
axial U(l) anomaly. This anomaly in QCD prevents the 7/ meson from being a Goldston 
boson which is manifested in its large mass (m^ = 958 MeV), a mass which is much larger 
than the masses of other members of pseudoscalar nonet (see. App. [A}. Hence the rj' 
meson is not included explicitly in the conventional SU(3) ChPT, albeit its effects are 
hidden in coupling constants. A recent extension of ChPT methods has provided tools 
which enable to include rj' in a consistent way and perform reliable calculations. One of 
these tools is a chiral unitary approach, which is based on the chiral perturbation theory 
and the unitarization using the Bethe-Salpeter equation. 

Concernig hadronic decays of the rj and rj' mesons into three pion systems, we can 
certify that 3tt can be in isospin state only if a two pion subsystem is in I = 1 state. 
In case of the 7r 7r 7r° system the two pion can have ^ = 0, 1, 2 but coupling with the 
remaining pion to 1^ = is only possible if ^ = 1. However, the (7r°7r°)/ = i dose not 
exist (a corresponding Clebsch-Gordan coefficient is equal to zero [8]) and as a conseqence 
the decay rj(rj') — > 7r 7r°7r° has to violate isospin. 

In the case of the rj(rj') — > 7r + 7r~7r decay taking into account the Clebsch-Gordan 
coefficients [8] one can write [9]: 



where: 




Thus the full wave function for the 3tt system reads: 

(3vr) /=0 = ^1 [|vr + )|vr )|vr-) - |vr°)|vr + )|vr-) - \*+)\ir-)\«*) + 

|O|7T + >|7r ) " + K°)|7r-)|vr + )] . (2.5) 

This wave function is antisymmetric against any exchange of pions: 7r° <-> 7r + , 7r - <-> 7r + 
and 7T° <-> 7T~. In particular, by applying charge conjugation we have: 

C(3vr) /= o = -(3tt) 7= o. (2.6) 

This is in contradiction with C = +1 for rj(rj'). Therefore the decay rj(rj') — ► tt + it~tt 
should violate C or I. 
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On the other hand there exist a G operator which is constructed from the C parity and 
isospin I2 operators in the following way: 

Ce i7Th . 



G 



(2.7) 



The eigenvalue of this operator is given by = (— 1) / Ac, thus Xq = —1 for pions and 
\g = +1 for the 77 mesons. Therefore the decay 77(77') — ► tttttt does not conserve G. 

Historically these decays were considered as electromagnetic processes with partial 
widths smaller than the second order electromagnetic decay 77(77') — ► 77. But it appears 
that the electromagnetic contribution is small [10,11] and instead the process is expected 
to be dominated by the isospin-violating term in the strong interaction. Due to this fact 
we can neglect the electromagnetic terms, and thus the decay amplitudes of the 77(77') — ► 37r 
decay becomes directly proportional to the quark mass difference (m^— m u ) [6]. The lowest 




Figure 2.1: Isospin violation in 77(7/) 
(courtesy of A. Kupsc [9]). 



37r decay. Lowest order effective lagrangian contribution 



order contribution to the decay mechanism is given by the Current Algebra (CA). Figure l2~Tl 
shows graphs consisting of a combination of the 7/ — 7r° mixing (a,b) and crossed graph (c) 
of the elementary low energy QCD processes - scattering of two pseudoscalar mesons. The 
partial width of the 7/ — ► tt + tt~tt° decay calculated using Current Algebra is 66 MeV [12], 
which is much below the experimental value of 294±16 MeV [8]. Gasser and Leutwyler had 
corrected this value to 160 eV [13] using the second order in the low energy expansion of 
the effective QCD Lagrangian. This change implies the importance of the 7T7t interaction 
in the final state and involving higher loop calculations should improve this value since 
they give a better description of the tttt final state interaction. An other approach which 
includes a pion-pion interaction up to higher orders uses the dispersion relation, which 
connects the imaginary part of the decay amplitude with the amplitude itself. There are 
two estimations based on this method [14, 15] but using different formalisms. They lead 
consistently to an enhancement of the decay rate by about 14%. 

It was suggested by H. Leutwyler that the decay width of the 7/(77') ~ 
sensitive to the mass difference of the light quarks following the relation [6]: 

r^/^+^o oc T (m d 



7T + 7T 7T° is 



777,, 



(2.8) 



More explicitly the decay width can be written in a convenient form as: 

1 



■niv')- 







Qd 
Q 



(2.9) 
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and the u and d quark mass 



where the dependence of the d and u quarks mass difference is contained in the Q term: 



1 

~Q 2 



m:, 



(2.10) 



and Qo is given by a relation: 



Q 2 



m 



D 



K 



m 



K 



ml 



ii i 
m„j.m i — m n 

K+ 7T+ 7T U 



(2.11) 



Qd and To correspond to Q and T calculated in the Dashen limit [16] where the quark 
masses are constrained by the assumption that the electromagnetic mass difference for 
kaons and pions are equal. Using the leading order expression of masses of the pseudoscalar 
mesons and applying the Dashem theorem, it was numerically calculated that Qd = 24.1. 
The decay width T is very sensitive to the exact value of Q and thus provides the precise 
constraints for the light quark mass ratios. Q determines the major axis of the ellipse of 
the light quark mass ratios which is given by the formula: 



m u \ | 1 / rn s \ _ 
m d ) Q 2 \m d 



The ellipse is shown in Fig. 12.2 



(2.12) 








0.5 



l m<i 



Figure 2.2: Ellipse which characterises the quark mass ratios (for the description see text). The 
figure is adapted from [6]. 



For the determination of Q we can use T$ together with the experimantal value of T. 
Recently Q value of 22.8 was derived by means of the dispersion relation approach from 
preliminary KLOE data on the rj — ► 7r + 7r~7r° decay [17]. In Fig. 12.21 the quark mass ratio 
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is indicated by the shaded uncertainity ellipse. The upper and lower sides of the ellipse 
is bound by two dashed lines corresponding to Q=22.7 ± 0.8. Prom the left it touches 
the hatched region, excluded by the bound Am > 0, calculated in the leading term in the 
expansion of powers 1/N C and extended to first non-leading order. Other bounds are given 
by the limits of the branchning ratios R = T^i^^o /T^i^^. The dot corresponds to 
Weinberg's value, and the cross represents the estimations described in the reference [18]. 
The Weinberg ratio corresponds to Am = 0, and is located at the boundary of this region. 
In particular this implies m^/m^ > \ and excludes a massless u quark. 

The Tq factor can be calculated in the chiral limit where the mass of quarks tends 
toward zero m u = = 0. Estimations of the decay width by the Chiral Perturbation 
Theory are based on the leading term of the expansion in the quark masses and the precise 
calculations of To in the isospin limit. In case of this method the normalization for the 
decay width must be obtained from other experminents with the electroproduction of the 
7/ meson. 

Gross, Treiman and Wilczek [5] have proposed another method of establishing the 
light quark mass difference by finding the ratio of partial widths for the isospin violating 
7/ — > ir + TT~ir° decay to the isospin conserving process of rj — ► i]tt + tc~: 

— >7r+7r-7r° , s 3 / 7TI,] — 777-,, \ 2 

-- v ~ (16.8)— ( — , (2.13) 



T, n ,^ m+n - 16 V rn^ 

where the factor 16.8 denotes the ratio of phase-space volumes. Formula (|2.13p was claimed 
to be true under two assumptions: that the amplitudes for both decays are constant over 
the phase space, and the amplitudes are related via [5]: 

A(t]' -> 7r + 7r~7r°) = sine A{rf -> r/vr+vr"), (2.14) 

where sin indicates the 7r° — rj mixing angle: 

S [ n e = - A — 2 «. 2.15 

4 m s 

The second assumption directly implies that the decay rf — ► 7r + 7r~7r° proceeds entirely 
through the channel rj' — ► 777r + 7r~ followed by 7r° — 7/ mixing. 

From equation (|2.13p one can see that the measurement of the ratio does not require 
the information from other experiments for normalization of the partial decay width. Addi- 
tionally with the simultaneous measurement of these decays, due to the similar final state, 
many systematical uncertainties will cancel. 

But recently Borasoy et al. [19] claimed that the light quark masses can not be extracted 
from the ratio (12. 13ft . They studied the two assumptions which were mentioned above using 
the U(3) chiral unitary framework which is in good agreement with rj' data as regarding 
widths and spectral shapes. They showed that results from the chiral unitary approach are 
in disagreement with these two assumptions. Concluding that on the theoretical side a still 
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and the u and d quark mass 



more sophisticated treatment including the final state interaction is required for univocal 
statements about the quark masses [20]. 

In our experiment we intend to establish the ratio r„/ _ +7r + 7r - 7r o/r„/_ >7r + 7r -„ hoping that 
the near future progress on the theory side will permit us to determine univocally from 
this ratio the light quark mass difference — m u . Further more, in combination with 
experiments determining values of quark mass ratios m s /md and m^/m^, it would be 
possible to obtain the absolut values of the u and d quark masses. 



3. Measurement method of the T^^+^-^o 

The partial width T r j_ yK +^- 1! o may be determined from the branching ratio of the r\ — > 
7T + 7r~7r decay and the total width of the if meson. In the WASA-at-COSY experiment 
the rj' meson will be created in two proton collisions via the pp — > pprj reaction. The 
measurement have to be performed at best possible conditions, because the rj — > 7r + 7r~7r° 
decay channel is very rare. So far the rj' — > 7r + 7r _ 7r decay was never observed and only 
an upper limit of 5% for the branching ratio has been established [4]. Recent theoretical 
calculations, based on a chiral unitary approach, predict the branching ratio of about 
1% [19]. 

3.1 Branching ratio BR(rj' — > n + n^iv ) 



The branching ratio is a quantity which informs about the probability of the particle 
decay into a specific channel. It is defined as: 

BRi = (3.1) 

J- tot 

where Tj and Y to t denote the partial and total (natural) width of the particle, respectively. 
Index i indicates a decay channel. From (13. ip we can see that the determination of partial 
widths will rely on the precise determination of branching ratio and of the total width. 
For the rj — > 7r + 7r~7r° decay the branching ratio can be explicitly expressed as: 

BR(ri'^^-^) = T ^: t + o r° , (3.2) 

where r, J /_> 7r +, r - 7r o is the searched observable. The T^j is known from the Particle Data 
Group (PDG) estimations to be T l °j = 0.202 ± 0.016 MeV/c 2 [80, and BR(?/ -» vr+vr-vr ) 
still remains to be established experimentally. With the WASA-at-COSY facility we plan 
to determine this branching ratio as the ratio of the number of events where rj' decays into 
the ir + 7r~ir° system (iV 7? /^ 7r + 7r - 7r o) to the number of all produced rf mesons (iV*?*): 

*- + "») = %^. (3.3) 

v' 

For this purpose the rj' mesons will be produced in collisions of the proton beam with 
the hydrogen pellet target via the pp — > pprj' reaction. The outgoing protons and the 

1 The COSY- 11 collaboration is working on a more precise estimation of this value with a dedicated 
measurement of the rj' mass distributions by using the missing mass techniques applied to the pp — > ppX 
reaction very close to the kinematical threshold [21]. 
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products of the decay of the rj meson will be detected and identified using the WASA-at- 
COSY detection system. More details of the measurement techniques will be given in the 
following sections. 

3.2 Description of the WASA-at-COSY detector facility 

The COoler SYnchrotron (COSY) is a storage ring operating at the Research Cen- 
tre Julich (Germany) since 1993 [22]. It delivers unpolarized and polarized proton and 
deuteron [23] beams in the momentum range between 300 and 3700 MeV/c. The first step 
of the particle acceleration takes place in the isochronous cyclotron (JULIC). Next the 
beam is injected into the 184 m long COSY ring (see Fig. 13. ip where the particles experi- 
ence further acceleration and at present may be used for experiments with the WAS A [1], 
ANKE [24] and TOF (external beam) [25] detection setups where the energy of the beam 
allows for production of all basic pseudoscalar and vector mesons. 




Figure 3.1: Schematic view of the COoler SYnchrotron ring in the Research Centre Julich. 
The presently used detector systems: WASA [1] and ANKE [24] and completed experiments 
COSY-11 [26] and PISA [27] are shown. 

The COSY accelerator is equipped with two types of beam cooling systems: an electron 
and stochastic cooling used for low and high energies, respectively [28]. Both cooling 
systems allow to decrease the momentum and geometrical spread of the beam. The whole 
process for the beam preparation from injection till final state of acceleration takes a few 
seconds. The ring can be filled with up to 10 11 particles, and the life time of the circulating 
beam ranges from minutes to hours depending on the thickness of the used target. 

In 2006 the WASA detector [29] has been transferred from the CELSIUS [30] facility in 
Uppsala to the COSY ring [22] in Julich where it was successfully installed and brought into 
operation. The WASA-at-COSY, shown schematically in Fig. 13.21 is a large acceptance 
detector which consists of three main parts: the Central Detector (CD), the Forward 
Detector (FD) and the pellet target system [31-33]. 

The Central Detector (CD) [34] is used for detection and identification of charged and 
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Figure 3.2: Cross section of the WASA-at-COSY detector with new FRH4 and FRH5 layers, 
mounted in August 2006 [43] and newly assambled two layers of the FWC [37]. Names of the 
detectors are explained in the text. 



neutral particles (7, ir + , tt~ , e + , e~) which are the decay products of short lived mesons like 
7T°, r] and if. It covers scattering angles between 20° and 169°. The most outer part of the 
CD constitutes an electromagnetic calorimeter (SEC) built out of 1012 Csl scintillating 
crystals positioned around the interaction point [35]. The crystals are ordered in 24 rings 
between the iron yoke and super-conducting solenoid which provides an axial magnetic 
field enabling the determination of the momenta of charged particles by measuring their 
tracks with the Mini Drift Chamber (MDC) [36]. The MDC is mounted around the beam 
pipe inside of the solenoid. It is a cylinder which consists of 17 layers of straw tubes. Each 
straw is made out of thin mylar foil tube, with a gold plated sensing wire in its center. 
The drift chamber is surrounded with the Plastic Scintillator Barrel (PSB) [36] used for 
trigger purposes and for the determination of the energy loss for charged particles. PSB 
together with SEC and MDC permits the identification of charged particles by means of 
the energy loss method (AE). 

The forward part of the WASA detector is built out of thirteen scintillating layers and 
four layers of straw drift chambers. The forward detector was designed for the detection 
and identification of protons, deutrons and He nuclei. The detector closest to the scattering 
chamber - Forward Window Counter (FWC) - covers the conical exit window of the axially 
symmetric scattering chamber towards the Forward Detector assembly [37]. Presently it 
is used for trigger purposes, but in the future it will also serve as a start detector for the 
time-of-flight determination. 

For the track reconstruction the Forward Proportional Chamber (FPC) [38] is used. It 
provides information about particle scattering angles with a precision better than 0.2° [39]. 
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The FPC planes are rotated by 45° with respect to each other. Directly behind the straw 
detector the three layer Forward Trigger Hodoscope (FTH) [40,41] is placed. It consists 
out of 96 individual plastic scintillator elements arranged in three layers: two layers with 
modules in the form of Archimedean spiral and one layer with cake-piece shaped modules. 
The FTH provides information about charged particle multiplicities in the FD used for 
the first level trigger, as well as in the offline track reconstruction [42]. Future trigger 
developments aim at a complete real time scattering angle reconstruction of individual 
tracks. This information, combined with the information of deposited energy in successive 
detector layers, will allow to determine the missing mass of forward going particles on the 
trigger level, and thus will provide a very efficient meson tagging [42]. The FTH detector 
can also be used to deliver the time information and it is thus helpful for the application 
of the TOF technique. 

The next five thick planes are called Forword Range Hodoscope (FRH) and are made 
out of cake-piece shaped plastic scintillator modules (see Fig. 13 .3|) . This detector enables 
identification of charged particles from the absolute value and the pattern of energy de- 
posited in different layers. Among the second and third layer of the FRH the Forward 



Figure 3.3: View of the Forward Range Hodoscope build out of the cake-pieced scintillaiting 
modules arranged in five layers [43]. 

Range Interleaving hodoscope (FRI) [44] is mounted. It is composed of two interleaving 
planes of 32 plastic scintillators bars aligned horizontally and vertically. The main function 
of the FRI is the measurement of the hit position of the charged particles improving the 
resolution for the vertex reconstruction. In addition, the FRI delivers also time information 
which can be used in the future application of the time-of-flight method. 

The Forward Veto Hodoscope (FVH) [45] is the last detection layer. It consists of 12 
horizontally placed plastic scintillator modules equipped with the photomultiplayers on 
both sides. The FVH gives the position of particle hits, reconstructed from time signals 
determined on two side of the module. It is important to note that this detector can also be 
used to deliver time information (t s t op ) for the TOF technique. It is also worth mentioning 
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that in the near future it is planned to extend the FVH by a second layer of vertically 
arranged scintillators [46]. 

Finally, depending on the studied reaction, a passive absorber (FRA) made of iron can 
be positioned in front of the FVH, enabling to disantangle between slow and fast particles. 
The thickness of the material can be choosen from 5 mm up to 100 mm. 

From the beginning of its operation at CELSIUS and COSY the WASA detector is 
equipped with a pellet target system [31-33], providing a stream of frozen hydrogen droplets 
with a diameter of about 35 /im. Pellets are passing to the interaction region through a 
thin 2 m long pipe. Usage of this type of target enables to achieve densities of up to 10 15 
atoms/cm 2 resulting in luminosities of up to 10 32 cm _2 s _1 when combined with the COSY 
beam of 10 11 particles stored in the ring. 

3.3 Identification of the pp — > pprj' — > pp^n^n — > ppir + n~^ reaction 
chain 

The proposed experiment aiming at the determination of the branching ratio for the 
rj — > tt + tt~tt° decay will be based on the production of the 7/ meson in proton-proton 
collisions and registration of its decay products. The complete reaction chain which needs 
to be identified (production and decay) reads: 

pp —> ppt] 1 —> ppTr + Tr~~ir — > ppir^ 'ir~ '77 . 

The measurement relay on the registration of all particles in the final state and on the 
determination of their four-momenta: 

In order to measure the energy of the two forward scattered protons the FRH1-5 planes 
in the Forward Detector are used which provide information of the energy losses. To 
determine the direction f of these two protons we use the straw chambers of the FPC and 
the FRI detector. Charged pions scattered under angles larger then 18° will be registered 
in the Central Detector. The MDC will provide information about the direction r(9, (ft) 
and the SEC will allow to measure the energy losses. In the case when the charged pion 
will be scattered forward (angular range of 2.5° - 18°) it will be detected by the Forward 
Detector in the same way as the protons. The two gamma quanta originating from the 7r° 
decay will be recorded in the SEC from which we will receive the information about their 
energy and direction. 

3.4 Invariant and missing mass techniques 

At the first stage of the data analysis we have to select from the full data sample 
only these events where tt + tt~tt° were produced and in the second step using a missing 
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mass technique we will determine fractions of selected events corresponding to the direct 
production and to the decays of rj meson. 

In order to identify signals from the 7r° mesons, we will reconstruct the invariant mass 
of two gamma quanta: 



which should be equal to the mass of the neutral pion within the expected resolution. Next 
we identify the charged pions on the basis of the energy measured in the SEC (E n ±) and 
their momenta (p w ±) derived from the curvature of the tracks reconstructed from signals 
measured with the MDC. In the case of the tt^ 1 the energy and momenta should fulfill 
(within the expected resolution) the relation: 



where m^± denotes the mass of a charged pion. 

Further on, from measuring the energy loss in the five scintillator layers of the FD 
and tracks in the FPC (in the future also from the time-of-flight) we reconstruct the four- 
momentum vector of the forward emitted protons [47]. Next in order to check whether the 
identified pions originate from the decay of the rj meson, we calculate a missing mass of 
the pp — > ppX reaction according to the equation (see Appendix IB]) : 



where mx corresponds to the mass of an unobserved particle, Eb eam ,pb earn denote the 
energy and momentum of the beam respectively, and E pi , E P2 ,p±,p2 represent energies 
and momenta of two registered protons. At present the experiment is still in preparation 
therefore as an example of the reconstructed missing mass distribution in Fig. 13.41 we show 
data derived from the COSY-11 measurement carried out at an excess energy of Q=15.5 
MeV above the threshold for the rj meson production [48]. The number of events which 
correspond to the production of the rj meson is marked by the letter S, and the field of 
the slice under the rj peak correspond to the number of the background events marked 
by B. The background originates from the direct two, three and more pions production. 
It maybe treated as the upper limit for the expected background of the pp — > ppir + iT~Tr 
reaction. In the WASA-at-COSY experiment we expect a similar spectrum for the missing 
mass of the pp — > pp(X = tt + tt^tt°) events. However, the signal to background ratio 
could be different due to the different mass resolution of the WASA-at-COSY detector. 
Additionally it will vary with the excess energJE Thereforee, in order to find an optimum 
beam momentum for the BR(r/' — ► 7r + 7r~7r°) determination the energy dependence (i) of 
the signal (pp — ► pprj — > pptt + tt^tt°), (ii) of the background (pp — ► pp7r + 7r~7r ), (hi) of 
the missing mass resolution and (iv) of the detection efficiency has to be established. The 
upper limit of number of background events may be obtained from results achieved by the 

2 Excess energy Q is a kinetic energy available in the reaction exit channel in the center-of-mass system. 




(3.4) 




(3.5) 




(3.6) 
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Figure 3.4: Example of the missing mass distribution for the reaction pp — ^ ppX from the 
COSY-11 measurements at the excess energy Q = 15.5 MeV [48]. S denotes the signal originating 
from the rf production and B indicates the background under the peak. 



COSY-11 experiment (where one spectrum is shown in Fig. I3.4[) . A detailed description of 
the derivation of the dependence of S(Q) and B(Q) will be given in the next chapter. 
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4. Conditions for the determination of the BR(r/ — > 7t + tt tt°) 
via the pp —> ppr[ —> pp7T + 7T^7T° reaction chain 



The number for all produced r/ mesons iV*°' (see eq. 13. 3|) can be determined from the 
known total cross section for the pp — > pprj' reaction [49-54] and the luminosity which will 
be established by the measurement of the reaction with the well known cross section e.g. 
from elastic scattering of protons (pp — > pp) [55]. 

When we assume that the relative statistical error of N^° l can be neglected due to large 
number of measured pp — > pp or pp — ► pprf events the formula for the relative error of the 
branching ratio is reduced to: 



a {BR) <t(N, 



rO 



^IZUE^ (4.1) 



BR iV r? /_ >7r + 7r - 7r o 

where cr(iV r; /_ >7r + 7r -, r o) denotes the statistical uncertainty of the signal. Assuming that the 
shape of the background is known the statistical error of the signal can be approximated 

as: 



ffOfy-Mr+ir-wo) ~ y N v' *° + N B, ( 4 -2) 

where Nb indicates the number of all background events under the signal. Equation 14.21 
was obtained under the assumption that the statistical errors propagate according to the 
formula [56]: 



a(Y(xt, ...,x n )) 



( dY , 



where the variables independent and a(xi) denotes their uncertainty. Thus 

the relative accuracy of the branching ratio can be expressed as: 

a(BR) = VAy^+^o + N B 
BR iV r? /_ >7r + 7r - 7r o 

The number of events where r( decayed into three pions which we expect to register with 
the WASA-at-COSY facility, can be expressed by the formula: 

At 

N S (Q) = Nj^n-AQ) = offiQ) ■ B R(V' ^ ^vr-vr ) • A(Q) f L ■ dt, (4.5) 

o 

where L is the luminosity which is expected to be around 10 32 cm _2 s _1 , cr*?* denotes the 
total cross section for the rj meson production which depends on the excess energy Q, 
BR is the supposed branching ratio for that decay, At indicates the measurement time 
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and A(Q) representes the acceptance of the WASA-at-COSY detector for the measured 
reaction. 

In the following section we will estimate the uncertainty of the branching ratio deter- 
mination as a function of the excess energy and measurement time taking into account the 
energy resolution of the forward part of the WASA detector. To this end, in the following 
we will parametrize the energy dependence of the total cross section for the rf meson and 
for the multimeson production. Presented estimations of the multimeson production will 
allow to estimate an upper limit of the background for the rf — > -k + tt~-k° decay. This is 
because it will be based on the missing mass distributions which includes also other pion 
channels like e.g. pp — > pp2ir, pp — > ppAir, pp — » pp2iri]. 

4.1 Parametrization of the total cross section for the rf meson 
production 

In order to estimate the expected production rate of the rf meson given by formula 
(|4.5p it is mandatory to know the total cross section for the production of this meson. 
From measurements of the COSY-11, DISTO and SPESIII [50-54] collaborations we have 
several experimental data points for the excess energy range from 1.3 MeV up to 150 MeV. 
The data are shown in Fig. 14.11 where it is clearly seen that near the threshold production 
of the rf meson depends strongly on the excess energy and hence this dependence must be 
taken into account. 



S 




10 

Q[MeV\ 

Figure 4.1: Experimental data points of total cross section for the pp — > ppvf reaction de- 
rived from the COSY-11, DISTO and SPESIII measurements [50-54]. The solid line denotes 
parametrization of the cross section using formula l|4.7p . The figure was adapted form refer- 
ences [48,60]. 



Generally the production cross section can be described as an integral of the square of 

ps 



the production amplitude |M| 2 over the phase space volume V ps using formula: 



off(Q) = jj dV ps \M\ 2 , (4.6) 

where F denotes the flux factor of the colliding particles. This energy dependence of 
the cross section can be written in a closed analitical form using the Faldt and Wilkin 
model [57,58] which takes into account the proton-proton Final State Interaction (FSI): 

= j= = C 2 Q2 2 * (4.7) 

where C\ and C2 denotes the normalization constant, e stands for the binding energy 
[57,58], and the A(s,mf,m|) is the triangle function [59], defined as: 

X(x, y, z) = x 2 + y 2 + z 2 - 2xy - 2xz - 2yz. (4.8) 

The free parameters e and C2 have to be established by fitting the function (|4.7|l to the ex- 
perimental data. The estimation of these parameters and the fit was done in reference [60], 
and the obtained values amount to [48]: 

e = 0.62 ±0.13 MeV 

Co = 84 ± 14 mb. 



Knowing the parametrization from (14. 7ft and having the values of the parameters we can 
compute the cross section for the 7/ meson production for excess energies near threshold. 
The result of the fit together with the experimental data is shown in Fig. 14.11 

4.2 Differential cross section -j 2 - 1 for multimeson production 

dm \ m=m rj i 

For the estimation of the background from the direct 7r + 7r~7r° production in the miss- 
ing mass spectrum of the pp — > pp(X = tt + tt~7t°) reaction, expected to be observed by 
the WASA-at-COSY detector, we need to calculate the da (pp^PP^ w n ) anc j the 

x m x =m ri i 

missing mass resolution (A w ) of the WASA detector setup. 

To our knowledge there are no data available on the invariant mass distributions or even 
on the total cross section for the 7r + 7r _ 7r° production in the proton-proton reaction near the 
kinematical threshold of the rf meson production^. Therefore, in order to estimate at least 
an upper limit of the differential cross section for the 7r + 7r~7r° production we have taken the 
missing mass spectra of the pp — > ppX reaction determined by the COSY-11 collaboration 
for several beam energies near threshold for the rj meson production [49-52] (see e.g. Fig. 



1 The total cross section of the pp — > ppn + Tv tt° reaction has been established experimentally only at 
three proton beam energis [61-63] (see App. [D[) . 
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13.41) . The differential cross section for the background was calculated according to the 
formula: 

da B _ N B (Q)(rf(Q) 



dm N S (Q) Am ' 
which was obtained by dividing Ng(Q) and Nb(Q) as given below. 
The number of all measured pp — > ppr/' events can be expressed as: 

At 



(4.9) 



N S (Q) = of(Q) ■ A(Q) J L-dt, 



(4.10) 



and the number of the measured background events can be approximated by: 

At 

A(Q) J L-dt, 



N B (Q) ~ ^ • 



dm 



(4.11) 



where A denotes the acceptance of the COSY-11 detector, which in a very good approxi- 
mation depends only on the mass of the produced system and on the excess energy Q [64], 
L indicates the luminosity, At is the time of the measurement and Am denotes the range 
of the missing mass values around the rj' signal which is approximated by: 



Am w 2F w 4.7a std . 



Figure |4~2] and Tab. HJ] show values of the 
data [50-52] by means of formula l4~9l 



(4.12) 

extracted from the experimental 




10 
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Figure 4.2: Inclusive differential cross section for multimeson production derived from the 
COSY-11 data [50-52]. The superimposed line shows the function = a ■ Q 13 fitted to the 
COSY-11 data, treating a and (3 as free parameters. 
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Table 4.1: Differential cross section -j 3 — for multimeson production in proton-proton 

collisions extracted from the COSY-11 data [50-52]. The systematical error amounts to 15% as 
established in reference [52]. 

A good description of the data was obtained by the function of the form: 

da 



dm x 



(Q) = a-QP, (4.13) 



where a and [5 are free parameters which, for Q expressed in units of MeV, were estimated 
to be a = 0.64 ± 0.14 nb/MeV and j3 = 1.662 ± 0.081 [7,73]. Therefore, conservatively 
the signal to background ratio for the WASA-at-COSY detector can be calculated using 
formula [4791 by replacing Am by the WAS A missing mass resolution, and by replacing the 
cross section cr*?* by the product a*?* x BR(rf — > 7r + 7r~7r ). 

In order to understand the production mechanism of the 7r + 7r _ 7r° system more precise 
studies are needed, because till now it did not receive a proper attention neither experi- 
mentally nor theoretically [7]. Therefore, we have conducted investigations by simulating 
the distributions for several production mechanisms and compared the results to the 
values of -jjjz- extracted from the experimental data. A more detailed description 

x m x =m v , 

of these studies is presented in appendix [Fj 



4.3 Energy resolution of the Forward Detector 

The last parameter needed to be calculated for the estimation of the background 
expected to be observed with WASA-at-COSY is the missing mass resolution Am w which 
depends on the energy resolution of the Forward Detector. 
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The Forward Detector was previously used for studies of the protons from the pp — > pp?7 
reaction. It consists of plastic scintillator layers which measure the energy losses, and on 
their basis the kinetic energy of the particles can be reconstructed. The thickness of the 
detector was optimized for a measurement of protons with kinetic energies in the range 
from 100 to 550 MeV. 
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Figure 4.3: Relative energy resolution as a function of the proton kinetic energy T achieved using 
different numbers of FRH planes of the WASA-at-COSY detection setup (courtesy of H. Calen [66]). 

In order to study the production and decays of the 7/ meson an upgrade of the Forward 
Range Hodoscope was necessary, due to the higher background to signal ratio and higher 
energies of outgoing protons from the pp — > pprf reaction (kinetic energy ranges from 300 to 
800 MeV). The extension was done by removing one 11 cm thick layer and instead adding 
two new layers of 15 cm thickness each built out of cake-piece shaped plastic scintillators 
with a photomultiplier tube attached at the end [43]. At present the total thickness of 
FRH detector amounts to 63 cm, which enables to stop protons with kinetic energies up to 
360 MeV. Thanks to the new layers the accuracy of the determination of proton energies 
was improved by about 25%. Fig. 14.31 shows the relative energy resolution of protons for 
different number of FRH planes. The accuracy of the energy reconstruction is a crucial 
point in the analysis. The systematical and statistical errors in the evaluations of branching 
ratios for the studied decay channels will strongly depend on the precision of the missing 
mass reconstruction, which in turn will depend on the precision of the reconstruction of 
the momenta of the forward emitted protons. The effective energy resolution of the newly 
assembled Forward Range Hodoscope is estimated to be about a = 3% for protons from 
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Figure 4.4: (left) Accuracy of the missing mass reconstruction as a function of fractional kinetic 
energy resolution simulated for a beam momentum of Pb eam =3.35 GeV/c. The ^ denotes the 
beam momentum resolution and the 'Std vtx' indicates the standart vertex resolution, (right) 
Simulation assuming = 0.1% and standard vertex resolution, for three beam momenta as 
indicated in the figure (courtesy of A. Kupsc [9]). 

the pp — ► pprj' reaction, at a beam momentum of 3.35 GeV/c. 

In order to investigate the missing mass resolution for the pp — ► ppX reaction we 
assumed that the COSY beam momentum spread is Ap/p « 10~ 3 , perpendicular beam 
profiles: horizontaly <tx=2 mm, vertical cry=5 mm [67] and the pellets are passing through 
the interaction point distributed homogeneously in a cylinder with a diameter of 2.5 mm. 
The result is presented in Fig. 14.41 where the width of the missing mass peak is plotted as 
a function of the relative resolution of the kinetic energy. It is seen that even in the case 
of a perfect energy resolution of the FRH the missing mass resolution amounts to about 
4 MeV/c 2 . If we further assume a perfect interaction point (vertex) determination the 
contribution from the beam momentum spread is seen to be at the order of 1 MeV/c 2 . In 
the right panel of Fig. 14.41 the dependencies were plotted for three different beam energies: 
3.24, 3.35 and 3.5 GeV/c from the near threshold rj production region. The broadening 
of the signal with increasing beam momentum is a kinematical effect due to the error 
propagation discussed in detail e.g. in reference [68]. 

Taking into account all effects which we mentioned above, the missing mass resolution 
can be parametrized as a function of excess energy and the kinetic energy resolution of 
forward detector by the following formula [69]: 

7 ATI 

T = (0.87 + 1.25)(0.48 + —0.17). (4.14) 
This relation will be useful in the next section where the accuracy of the branching ratio 
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Figure 4.5: Acceptance of the WASA-at-COSY detector as a function of the excess energy: 
(left) for the pp — > pp?7' — > pp7r + 7r~7r° and (right) for the — > pp7r + 7r~7r° reaction. Solid line 
indicates acceptance calculated assuming a homogeneous phase space distribution, and the dashed 
line denotes the acceptance where additionally a proton-proton FSI was taken into account. 



determination as a function of excess energy will be calculated. 

Additionally we have included the acceptance of the WASA-at-COSY apparatus for 
detecting the pp — > pprj' — ► ppir + ir~ir and pp — > ppir + ir~ir° channels. We have taken 
into account the geometrical acceptance of the central and the forward detector which for 
detecting the protons and pions covers the following ranges of the polar angle: 2.5° ^ 9 ^ 
18°, and 20° ^ ^ 169°. Furthermore the Final State Interaction (FSI) between outgoing 
protons was included using the square of the on-shell proton-proton scattering amplitude 
calculated according to the Cini-Fubini-Stanghellini formula including the Wong-Noyes 
Coulumb corrections [70-72] (for details see appendix [E]). Fig. 14.51 indicates the excess 
energy dependence of the acceptance assuming that the phase space is homogeneously 
populated (solid line) and including the FSI between protons (dashed line). 



4.4 Accuracy of the branching ratio determination 

The direct three pion production and the rf — > tt + tt~tt° decay will be disentangled by 
using the missing mass of the two outgoing protons measured in the forward detector. It is 
worth noting that in order to distinguish (on the base of a statistically significant sample 
of events) the direct and resonant multi-pion production it is mandatory first to select a 
sample of events with the studied final state channel (e.g. 7r + 7r~7r ) and then only for 
this selected sample to construct a distribution of the missing mass to the proton-proton 
system. 

To estimate the accuracy of the BR(7/ — > 7r + 7r~7r ) determination we have parametrized 
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the total cross section for the rf meson production (eq. l4.7p and have etablished a parametriza- 
tion for an upper limit of the background production (eq. 14.131) . Further on, the missing 
mass resolution taking into account effects which are related to the energy resolution of 
the forward detector and the beam and target spread (eq. I4.14p has been parametrized. 
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Figure 4.6: The relative accuracy of the determination of the BR(?/ — > 7r + 7r~7r°) as a function 
of the excess energy Q for the pp — > pp?]' reaction. One week of data taking with WASA-at-COSY 
with luminosity L = 10 32 cm~ 2 s _1 is assumed for the calculations. 

Applying these parametrizations and assuming five values of the BR(?/ — ► 7r + 7r~7r°): 
5%, 2%, 1%, 0.75% and 0.5% we have made calculations for the relative accuracy of the 
branching ratio determination. For the calculations we considered the range of values 
from the established upper limit of 5% down to 0.5%. Assuming for example that the 
measurement will last one week with a luminosity of L=10 32 cm _2 s _1 , we obtained the 
result shown in Fig. 14.61 where the relative error of the branching ratio is plotted as a 
function of the excess energy Q. We can see that the relative statistical error scales nearly 
with the value of the assumed branching ratio, and that the optimum accuracy is achieved 
for excess energies between 60 and 90 MeV independently of the BR magnitude. Relation 
(|4.4p implies that the relative statistical error of the branching ratio determination will 
improve with time as 1/Vi, as it is shown in Fig. 14.71 (left) for a beam momentum of 
Pbeam = 3.45 GeV/c corresponding to an excess energy of Q = 75 MeV. The plot shows 
that if the branching ratio was equal to 0.5% a relative accuracy of 10% would require two 
months of data taking [65]. 

Additionaly, for the case if the BR(?)' — ► 7r + 7r~r/) is to small to be observed we have 
estimated the upper limit at a confidence level of 90% as a function of measurement time. 
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Figure 4.7: (left) The relative accuracy of the BR(?/ — ► 7r + 7r~7r ) determination as a function 
of measurement time for three assumed values of the branching ratio and a beam momentum of 
Pb = 3.45 GeV/c. (right) Upper limit of determination of the BR(r/ — > 7r + 7r _ 7r°) for a confidence 
level of 90% versus the time of the experiment. Estimations are performed assuming a luminosity 
of L = 10 32 cm -2 s" 1 and beam momentum of Pb = 3.45 GeV/c. 

We can see that e.g. a sensitivity of 0.001 can be reached after few weeks of data taking 
with the WASA-at-COSY detector. As it was pointed out the theoretical calculation based 
on the chiral unitary approach predicts the branching ratio af about 1% [19]. From the 
right panel of Fig. 14. 71 we may infer that this prediction can be falsified or confirmed within 
a few days of beam time using the WASA-at-COSY detection setup. 

It is worth to mention that an additional source of background, not discussed here, 
comes form other decays of the rj meson involving similar particles: rj' — > 7r + 7r~7/ and 
r/ — ► uij. This background cannot be suppressed using the missing mass method, yet it 
can be identified using reconstructed invariant masses of the decay products. 



5. The Time-of-Flight method 



In order to study rare decays of the r\ and rf mesons produced in the reaction 
pp — > ppi](rj') using the WASA-at-COSY detection setup, one has to determine the four- 
momentum vectors of the decay products and of the forward scattered protons. At present 
the protons four-momentum vectors are derived only taking into account the energy loss 
measured in the Forward Detector. In this chapter we will present a possible future im- 
provement of the momentum reconstruction using the time information from the scintil- 
lators of the Forward Detector [74,75]. We will also describe the alghorithm which can 
be implemented to reconstruct the momentum of particles crossing the Forward Detector 
using time information obtained from the crossed detection layers. 

5.1 Proposal of using the Forward Detector for a TOF method 

In order to measure time-of-flight (TOF) of charged particles emitted in forward direc- 
tion we can use all scintillator detectors which are placed in the forward part of the WASA 
setup. But for the purpose of this thesis we will consider only the thin plastic scintillator 
detectors: Forward Window Counter, FWC (2x3 mm), Forward Trigger Hodoscope, FTH 
(3x5 mm), Forward Range Intermediate Hodoscope, FRI (2x5 mm), and the Forward 
Veto Hodoscope, FVH (20 mm). The five layer Forward Range Hodoscope can also be 
used to determine time, but it is made of scintillators which thickness ranges from 11 to 
15 cm, and more precise studies of the light signal generation and propagation would be 
necessary for the usage of the time information. 

Using only thin detectors, enables us to measure time in eight points on a distance of 
about 2.0 meters. But in practice one can combine the individual times measured in layers 
of each detector to obtain one time value per detector. Additionally the information about 
the trajectory of particles is needed to recontruct the time in the FWC, FTH and FRI 
counters, due to one side readout of the scintillating modules (only one photomultiplier) . 
The direction of particles will be obtained from measuring tracks in the straw chambers 
(FPC). The FPC chambers are placed between the first two scintillator detectors - FWC 
and FTH. The material which they are made of is very thin, therefore we can neglect the 
particle energy losses in these layers. 

The five layers of the Forward Range Hodoscope, FRH, are placed between FTH and 
FVH. Their total thickness amounts to 63 cm, therefore particles which travels through 
the FRH will be significantly slowed down. However, the deposited energy in the FRH 
layers is measured with good accuracy in the order of few per cent, and by combining the 
TOF information and the deposited energy one can reconstruct the initial velocity of the 
particle at the interaction point. Additionally the time signals will be measured by the 
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Figure 5.1: Schematic view of detectors used for future TOF determination with the distance to 
each other given in meters. The FVH2 is the second veto layer which is planned to be build in the 
near future [76]. 



FRI detector which is placed between the second and third layer of the range hodoscope. 

Generally, taking into account the time information from thin scintillators and the 
energy loss from thick scintillators we can write that all information which will be measured 
is collected in eight independent values of time and five energy losses: 

(t, AE) exp = t exp , t e 2 xp , t e 3 xp , t exp , t exp , t exp , t exp , t exp , AE exp 7 AE e 2 xp , AE e 3 xp , AE ex p , AE exp . 

(5.1) 

Using Monte-Carlo methods we can simulate the same quantities as from the experiment: 

I + A TT\ j.mc + mc + mc + mc + mc j.mc 4.mc ±vnc a rpmc \ Tpmc \ rrmc \ rpmc \ jpmc 

(E, l\Hj) mc =t 1 , t 2 , t 3 , t 4 , t 5 , t g , t 7 , Eg , , l\tL 2 1 5 Aii 4 J A - C/ 5 • 

(5.2) 

Therefore, in order to reconstruct the correct particle energy (which the particle had at 
the interaction point), we have to assume different true particle energies, and compare 
the Monte-Carlo predicted times and energy losses for each assumption with the measured 
qunatities until we find the optimum (most probable) value for the kinetic energy of the 
particle. A figure of merit in comparing this two samples could be a x 2 criterion. A more 
detailed description of the algorithm is given in section 15.61 

5.2 Description of the detectors used for the time measurement 

In this section we will briefly describe the detectors relevant for the time measurement. 
The first detector which delivers time information on the way of the forward scattered 
particles is the Forward Window Counter (FWC) [37]. The light collection in the detector 
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is optimized to keep the detection efficiency as homogeneous as possible over the full 
hodoscope area. The detector is 48-fold segmented and is composed of two layers a 24 
elements made out of 3 mm plastic scintillator. The first layer is of conical shape whereas 
the elements of the second layer are assembled in a plane. The elements of the second plane 
are rotated by one half of a module (7.5°) with respect to the first layer. This geometry 
provides a complete coverage of the forward area without holes. In addition the 48-fold 
granularity coincides with the 48-fold granularity of the FTH. We expect to achieve a time 




Figure 5.2: Schematic view (left) of the Forward Window Counter (FWC) [76] which is the 
"start" counter for TOF measurement (tijia), and (right) of the Forward Trigger Hodoscope 
(FTH) (t 3 ,U,t 5 ). 

resolution for each layer in the order of a (FWC) « 200 ps. A schematic view of FWC is 
shown in the left panel of Fig. I5.2L 

The next detector yielding time information on the way of the forward flying particles 
is the Forward Trigger Hodoscope (FTH). It consists of three scintillating layers which 
were renewed in the last year [42,77]. Setup has a highly homogeneous detection efficiency 
and shows a fairly uniform behavior [42]. In case of using the FTH detector to measure 
time it is expected to achieve a resolution of about cj(FTH) 200 ps for each layer. The 
front view of the FTH detector is shown in the right panel of Fig. 15.21 

FRI is a thin scintillator hodoscope, designed to provide fast spatial and time informa- 
tion from inbetween the FRH. It consists of two layers of horizontal and vertical scintillator 
bars of 5 mm thickness each, and a maximum length of 1405 mm. Each bar is read out on 
one side via small and fast photomultiplier with light guides of fishtail type. The expected 
time resolution of the FRI is about cr(FTH) 300 ps. Schematic view of the FRI is shown 
on left panel of Fig. 15.31 

The last thin detector in the pathway of a forward scattered particle is the Forward 
Veto Hodoscope (FVH). The FVH is a veto detector build out of 12 horizontal scintillator 
bars which are read out on both sides (right panel of Fig. 15. 3ft . The accuracy of a time 
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measurement amounts to about a « 150 ps. In future it is planned to build a second 
veto detector with 20 modules placed vertically and readout by photomultipliyers on both 
sides. This detector will be optimized for the TOF measurment, and it will be placed 
about 80 cm down stream with respect to the first FVH [46]. 




Figure 5.3: Schematic view of the: (left) Forward Interleaving Hodoscope(FRI) (^,£7), and 
(right) the Forward Veto Hodoscope (FVH) which will be the "stop" counter (ts,tg). 

By combining the time-of-flight method and presently used energy loss technique we 
expect to increase the momentum determination by a factor of about 1.5 and increase 
correspondingly the accuracy of particle identification. In order to estimate an expected 
improvement we have conducted Monte Carlo simulations of the pp — > ppr( reaction in 
the range of few tens of MeV above the threshold where the velocity of outgoing protons 
is around [5 « 0.75. The determination of the kinetic energy of protons from energy loss 
will at best be possible with a fractional resolution of ^P-=3%, which is equivalent to a 
relative momentum resolution of (see App. C): 

a M = l±HL^^ 2A% (5 .3) 
p T + 2m T y ' 

The determination of the time-of-flight in vacuum would lead to a fractional resolution of 
the momentum reconstruction expressed by the following formula: 

gfr) 1 °{tof) ,,,, 

p l-(3 2 tof ' 1 J 

Thus, from the above anticipated resolution of about 150 ps for FWC and FVH counters 
even with this two detectors, one could obtain a fractional momentum resolution of ^M- = 
3%. But this is only a conservative limit of the expected improvement since the particles 
will be slowed down in the FRH and the time will be measured by more detectors. 

5.3 Fractional energy resolution from the TOF measurement 



To check the posibility of energy reconstruction from the time-of-flight method we 
have performed studies of fractional energy resolution as a function of the kinetic energy 
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and kinetic energy as a function of time-of-fiight between start and stop detectors. We 
have simulated homogeneously in the phase space the pp — ► pprj' reaction using the GEN- 
BOD [78] procedure with a nominal beam momentum of Pbeam — 3.35 GeV/c without 
taking into account any beam spread. 
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Figure 5.4: The time-of-flight between the FWC and FVH detectors as a function of the proton 
kinetic energy. For the calculations the energy losses in the scintillators of the Forward Detector 
were taken into account. 

To calculuate the kinetic energy as a function of the time-of-flight we used only two 
detectors here: the FWC as start counter and the FVH as stop counter, which are about 
2 m apart (see Fig I5.ip . The energy loss in the detector material of FD was calculated 
using the Bethe-Bloch formula [79], which for plastic scintillators and protons in the kinetic 
energy range from 200 to 830 can be approximated by [8]: 



dE 5 
— — = const ■ p 3 , 
dx 



(5.5) 



where the factor const = 1.76 was calculated [69] from the experimental data [80]. 

The time-of-flight between the start and stop detectors was calculated iteratively using 
numerical Monte-Carlo techniques. Technically for each proton flying through the detector 
the time-of-flight was derived by summing very small time intervals for passing a distance 
on which the changes of the ionization power (dE/dx) can be neglected: 



TOF 



dx 

W) 



E(d) 



dE 



(5.6) 
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The plot in Fig. 15.41 shows the time-of-flight calculated as a function of the proton 
kinetic energy. We can see that with decreasing energy the time-of-flight between the start 
and stop couters increases. And for particles which have an energy smaller then 360 MeV 
the time-of-flight cannot be established because they do not reach the stop counter. 
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Figure 5.5: (left) Fractional kinetic energy resolution estimated for three considered techniques 
- energy loss, DIRC and TOF. (right) Distribution of kinetic energy of protons emitted from 
pp — ► ppr{ reaction simulated for a beam momentum of 3.35 GeV/c. 



Another possible improvement of the energy measurement of the forward emitted par- 
ticles assumes the installation of a DIRC detector which would enable to determine the 
velocity of protons [81,82]. Therefore, in the following we will compare the missing mass 
resolution obtained using the energy loss method with the resolution expected when us- 
ing the time-of-flight method and also with the resolution achievable when using a DIRC 
counter. 

For the comparison we expressed the properties of time-of-flight technique, DIRC and 
energy loss method in terms of fractional kinetic energy resolution. The result is shown in 
the left panel of Fig. 15.51 it indicates that both DIRC and time-of-flight can yield better 
resolution than that obtained from energy loss only. However, for the estimation of the 
missing mass resolution it must be taken into account that the TOF method can be used 
for protons which passed the whole detector with an energy T ^ 360 MeV, whereas the 
DIRC can deliver signals only for protons above the Cerenkov threshold which corresponds 
to 500 MeV [81]. This implies that only a fraction of protons (see right panel of Fig. 15.51) 
can be reconstructed by means of these two methods. The right panel of Fig. IBTBl illustrates 
that the energy of both protons can be reconstructed for only about 20% of events using 
the DIRC detector, and for about 70% of events using time-of-flight method. In other cases 
the energy of either one or both protons must be determineted using energy loss technique. 
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5.4 Accuracy of the missing mass reconstruction 

For a rough estimation of a missing mass accuracy we have assumed that we will 
measure the time only with the FWC and FVH but to account for the fact that in reality 
we will use much more detectors we assumed the precision of <r(tof) = 100 ps. Further 
more we have included the geometrical acceptance of the WASA Forward Detector for 
particles emitted in the forward direction, which is: 2.5° ^ 9 ^ 18°. The simulation did 
not include contributions due to a finite interaction region and the momentum spread of 
the COSY beam, since we focused on the comparison of the resolution resulting from the 
reconstruction of the proton four-momenta by the discussed methods. 

In order to calculate the missing mass distributions we have simulated proton four- 
momenta using an event generator based on GENBOD [78], which produces particle four- 
momenta in the CM-frame homogeneously distriuted in phase space. These particle four- 
momenta were then boosted into the lab-system by means of a Lorentz transformation. To 
reproduce the missing mass as measured in the experiment we then smeared the kinetic 
energy of the simulated protons applying the known energy resolutions for each of the 
studies reconstruction methods. Next using the calculated four-momenta of protons the 
mass of an unobserved particle was calculated according to the equation: 

m X = J {E beam + mtarget - E pi - E p2 ) 2 - {pbeam ~ Pi ~ P2) 2 , (5.7) 

In Fig. 15.61 the missing mass distributions for each method are shown. The left plot 
shows the result obtained using energy loss technique, middle plot corresponds to the DIRC 
technique and right to the time-of-flight method. 
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Figure 5.6: Missing mass distribution to the pp — » ppX reaction reconstructed under the as- 
sumption that the resolution of kinetic energy determination of protons will be as expected from 
energy loss measurement (left plot), DIRC (middle plot), TOF with a(TOF) = 100 ps (right plot). 

The simulations for the time-of-flight method and DIRC detector shows that the reso- 
lution of the missing mass reconstruction can be significantly improved in comparison to 
the results obtained when using only energy loss information. 
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5.5 Influence of the passive material on a precision of the energy 
determination 

New layers of the Forward Range Hodoscope (FRH 4 and 5) are placed in thin 8 mm 
plexi glass shield, to keep all scintillators modules in the correct position and to protect 
them against any mechanical damage. This material is not used in the measurement and 
we do not have information about the energy loss in this layers. We expect that in the 32 
mm (4 x 8mm) of the dead material the average energy loss will be about 10 MeV, and that 
the spread of this energy loss will be in order of 1 MeV or less. In order to estimate the scale 
of the effect we have conducted a simulations of a kinetic energy dependence of protons 
flying trought the Forward Detector as a function of the time-of-flight. In calculations 
we assumed a four different values for the error in calculating energy losses in the passive 
material: AE p i exi = 1 MeV, 2 MeV, 4 MeV, 10 MeV. The proton energy as a function of 
the time-of-flight for the assumed error is plotted in Fig. 15.71 where it is compared to the 
nominal dependence. We expect an error to be less then 1 MeV, however for the better 
visualization of the effect a much larger range was studied. 




TOF[ns] 

Figure 5.7: The kinetic energy of protons as a function of the time-of-flight. The lines indicates 
different assumption of the error for the calculations of the energy losses in the four passive layers 
of plexi. 

One can see from Fig. 15. 71 that the inaccuracy of the energy loss determination of about 
1 MeV will cause the error of kinetic energy determination in the order of fraction of MeV. 
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Therefore we can neglect this influence of the passive material on energy reconstruction. 



5.6 Computational algorithm for the reconstruction of particle 
momenta based on time signals in the Forward Detector 

In this section we will present the scheme of the algorithm developed for the recon- 
struction of particle energy which is based on time and energy losses measured in the 
Forward Detector. For the sake of simplicity further on we will consider only average times 
measured in each detector, and we will assume for simplicity that energy loss AE occurs 
in one solid block of scintillating material between the FTH and FRI and another block 
between FRI and FVH detectors, as shown in Fig. 15.81 The simplification is made having 
in mind that in general every detection plane can be treated as a separate detection unit, 
and that the derived equations can be easily generalized for more detection planes. 

In such case each event is characterized by a set of two values of energy loss measured 
in the FRH scintillation layers, namely: 



AE exp = (A£f p , AE e 2 x p ) 
and time information from four thin scintillators (FWC, FTH, FRI, FVH): 

,exp ,exp ,exp ,exp 



(5. 



(5.9) 




Figure 5.8: Schematic view of detectors used for the explanation of the algorithm to reconstruct 
protons energy based on time and energy loss information. 



The measured times can be expressed as a function of three parameters: 

t i = t i - 1 + f(E i _ 1 ,Al,AE), 



(5.10) 
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where E% denotes the particle energy after i-th detector, Al denotes the length of the 
particle trajectory between the two detectors, and AE indicates the energy loss between 
(i-l)-th and i-th time counter. 

Therefore we can calculate for each event a time using following equations: 



P(EoY 



5.11 



tr=tr+ l j^ y (5.i2) 

where /3(Eq) denotes the initial velocity of the particle. 

Having AE measured we can calculate the time in which particles pass a thick detector 
At(d) numerically. E.g. for the At(d,2) the formula reads: 

E 2 -AE(d 2 ) 

dE 



At(d) = / -= . (5.15) 

P(E)<§(E) 



E 2 



In order to reconstruct the true particle energy Eq which it had at the interaction point 
we compare simulated times i™ c with the measured times t^ xp . For the comparison we can 
construct the \ 2 according to the method of least squares: 



,,-exp tmc{Eo] 



X 2 (E ) = £ ( 3 ) (5.16) 



where, a\ indicates the resolution of the time mesurement in the i-th counter. One can, 
treat Eq as a free parameter, and determine the value of the particle energy when x 2 {Eo) 
reaches a minimum. To optimize the algorithm we may in practice scan the value of Eq 
in the range of uncertainty around the value established based on the energy loss method 
only, and using one of the minimalizing methods as e.g. the bisection method. 

The proposed method of particle identification using time signals, can be generalized 
for usage with real number of scintillator layers of the WASA-at-COSY detection setup. 



6. Summary and conclusions 



This thesis aimed at the estimation of the feasibility to study the branching ratio for 
the rf -> tt+tt-tt decay with the WASA-at-COSY detector. By measuring this isospin 
violating decay it is possible to derive the mass difference between the u and d quarks. 
This in combination with the ratios of the light quarks masses, can be used to derive the 
absolute value of quark masses, which cannot be observed directly. 

In order to established the optimum beam momentum for conducting the experiments 
with WASA-at-COSY setup we performed Monte-Carlo simulations taking into account the 
energy dependence of the signal and the background as well as the missing mass resolution 
and the detection efficiency. For the calculations we considered a range of values from the 
established upper limit of 5% down to a value by one order of magnitude lower (0.5%). As a 
result we found that the best accuracy for the measurement of the BR(?/ — > 7r + 7r~7r°) with 
the WASA-at-COSY detector is achieved when the rf meson is produced with an excess 
energy in the range between 60 and 90 MeV corresponding to beam momenta ranging from 
3.4 - 3.55 GeV/c. 

Further on we investigated the branching ratio uncertainty as a function of measurement 
time for one value of excess energy (Q = 75 MeV). As a result we can conclude that if the 
BR(r/ — » 7T + 7r _ 7r ) was equal to 1% as predicted based on the chiral unitary approuch [19] 
we would need five weeks of beamtime with a luminosity of 10 32 cm -2 s _1 in order to 
determine the branching ratio with an accuracy of 5%. We have established also the 
sensitivity for the estimation of the upper limit of the BR(?/ — ► 7r + 7r~7r°) in the case of 
no signal. We found that an upper limit on the branching ratio of 0.001 can be set on a 
confidence level of 90% after one month of beam time. 

Because there are no available data on the 7r + 7r _ 7r production in collisions of protons 
near the kinematical threshold for the rf meson production, we have used the COSY- 
11 data on the pp — > ppX reaction and established an upper limit for the background, 
expected to be observed with WASA-at-COSY detector. Using the missing mass spectra 
for the pp — > ppX reaction for several excess energies, the differential cross sections for 
the multimeson production was calculated. For further usage the excitation function was 
parametrized as a a • Q@ . By that means we can compute the values of the expected 
background as a function of the excess energy near the rf meson threshold. 

In the second part of the thesis we have proposed a possible improvement of the energy 
reconstruction of forward scattered protons using the time measured in the scintillator 
counters of the forward part of the WASA-at-COSY detector. 

The missing mass resolution based on a time-of-flight measurement was investigated 
and it was presented that it is similar to the resolution achievable with the energy loss 
method. Thus, combining this two independent techniques we expect to improve the 
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missing mass resolution by about a factor of 1.5. We also have presented the basic concept 
of an algorithm which will enable to determine the energy of particles based on the time- 
of-flight technique. In the near future we plan to implement the established algorithm in 
the WAS A analysis software. 



A. rj as a member of the SU(3) pseudoscalar meson nonet 



According to the quark model each (pseudoscalar) meson is a state of a quark and 
antiquark system. The three lightes quarks u, d and s, give nine possible qq combinations. 
In terms of the SU(3) symmetry we can write this as: 



3<8>3 = 801, 



(A.l) 



which includes an octet and a singlet. All possible meson states can be presented as a 
multiplet on a plot where the horizontal axis denotes the third component of the isospin - 
I3 and on the vertical axis the strangness - S is shown. Figure IA.1I presents a multiplet of 
the pseudoscalar mesons. 




7l u 




Figure A.l: Pseudoscalar meson nonet of the SU(3) symmetry group. 

In this picture we have pions: 7r + ,7r°,7r _ , kaons: K + ,K° and K~ ,K° and rj\ and rjg. 

The rj and 7/ meson are not pure states of the SU(3) pseudoscalar nonet but a combi- 
nation of a singlet 771 and octet r/g. The quark model give as wave functions for the pure 
states: 



V3 
1 



[uu + dd + ss) 



V6 



(uu + dd — 2ss) 



(A.2) 
(A.3) 



and from them the wave functions for physicaly observed states are constructed as: 

rj = rjg cos 9 — rji sin 9 (A. 4) 
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r\ = rjs cos 9 + rji sin 8 (A. 5) 

where 9 ps — 15.5°±1.3° [83] denotes the mixing angle which was established experimentally. 
Using equation (]A.4j) . (IA.5|) and the value of the mixing angle we obtain [84]: 

rj = 0.77 ^=(uu + dd)-0S3ss,. (A. 6) 

r{ = 0.63^(tm + dd) - 0.77ss (A.7) 
V2 

Both wave functions above are very similar and therefore the mass of these two mesons 
should in practice be also similar. However, the empirical mass of the rj meson is about 
two times larger than the mass of the 77 meson. 



B. The missing mass technique 



The missing mass technique is used to identify the particles which cannot be registered 
by the detectors due to short life time. We study the production of the rf meson via the 
pp — > ppX reaction where X denotes the unobserved particle. If we determine the four- 
momenta of all protons before and after the collision, the four-momentum conservation 
gives the relation: 

Pfceam + ^target = ^1 + ^2 + ^Xi (B-l) 

where Px corresponds to the four-momentum of the unobserved particle. 
The mass can be calculated from (jB.lj) as: 

m\ = \¥ x \ 2 = |P beam + F target - Fx - P 2 | 2 , (B.2) 

with P = (E,p) we can write: 

m X = {Ebeam + E tar get ~ E\ — E^) 2 — (Pbeam + Ptarget ~ Pi ~ P2f ' ■ (B.3) 

Using the formula above we can calculate the mass of the unobserved particle assuming 
that we can determine the four-momenta for all other particles participating in the reaction. 

In the laboratory where the proton beam collides with the hadrogen target the momen- 
tum of target is ptarget = an d energy is equal to the mass. Therefore in the laboratory 
system relation (1B.3|) can be written as: 

m X = ( E beam + m t arget ~ E X - E 2 f ~ (pbeam ~P\~ P2? ■ (B.4) 
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C. Relation between the fractional momentum resolution 
and the fractional time-of-flight accuracy 



To calculate the fractional momentum resolution related to the time-of-flight fractional 
resolution we use the formula for the momentum: 

p = m >yp p 2 = m 2 -f 2 p 2 , (C.l) 
where 7 denots the relativistic Lorentz factor 

7= ; 1 (C.2) 
and (3 indicates the velocity of the particle: 

P= 1 -, (C.3) 

with I being the distance for the time-of-flight measurement. Inserting the (IC.2j) and (IC.3|) 
to (jC.ip we get: 

P 2 = ^4- (C4) 



t 2 - I 2 " 

The standard deviation of the momentum distribution as a function of time resolution 
can be obtained by differentiating equation (|C.4j) : 

, 9n 2m 2 l 2 t . . s 

Applaying the relations t = jj, and cr(p 2 ) = 2pa(p) in eq. (IC.5j) we obtain a relation 
between the fractional momentum resolution and the fractional time-of-flight accuracy: 

-p--T^w~r- (C - 6) 

Further on, we will derive the relation between the fractional momentum resolution 
and the fractional energy resolution. By differentiating p 2 = E 2 — m 2 we obtain: 

2pa(p) = 2Ea(E), (C.7) 

and thus we can write: 



a(p) E 2 a{E) 1 a(E) 



(C. 



p p 2 E (3 2 E 

Prom formula (IC,8|) we can also estimate the relation of the fractional momentum 
resolution as a function of the fractional kinetic energy resolution. The total energy is 
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and the fractional time-of-flight accuracy 



given by E = T + m and by differentiating we can write: a(E) = <r(T), due to constant 
mass. Therefore we have: 

a(p) E 2 T a{T) 



p p 2 E T 
By further transformations we obtain: 

a{p)_ETa{T) ET a(T) _ ET a(T) 

p p 2 T E 2 — m 2 T (T + m) 2 — m 2 T 

a(p) ET a(T) E a{T) 



(C.9) 



(CIO) 



(C.ll) 



p T 2 + 2mT T T + 2m T 
Finally the fractional momentum resolution as a function of fractional kinetic energy res- 
olution reads: 

ff (p) = T + m a ( T ) (c l2 \ 

p T + 2m T \ ■ > 



D. Parametrization of the pp —> pp7T + 7T tt° total cross section 

The production process of three pions is not well known. Only a few experimental 
data points exist for the pp — > pp-Kinr total cross section and the dynamics of the process 
is still not well understood [7]. Therefore, for the estimation of the excitation function we 
have used a parametrization proposed by J. Bystricky [85] which is based on the expansion 
of the total cross section in a base of the generalized Laguerre polynomials. 
One can expressed the total cross section using the effective amplitude: 

Vp P ->ppn+7r-n° = \ F { X )\ 2 , (D.l) 

where the amplitude F(x) can be expanded into series of orthonormal functions L°(x): 

oo 

F(x) = J2 a nK(x), (D.2) 

n=0 

with 

L«(s)=e-isf£« (D.3) 
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Figure D.l: (left) The parametrization of the total cross section for pp — > ppir + ir~ir° and 
experimental data points from [61-63]. (right) Cross sections of the pp — ► ppij' reaction [51-54] 
compared to the total cross section for the pp — ► ppTr + ir~ir° process (solid line). 
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where C° denote the generealized Laguerre polynomials: 

C(x) = ±( n + a )t^l. (D.4) 

m=0 v 7 

We put 

x = c k \n(^-J (D.5) 

where is the beam kinetic energy at the reaction threshold, and c k is a constant, deter- 
mining the scale of T. 

One can see that the parametrization of the total cross section is in good agreement 
with the available experimental ponits for the three pion production. Considering now the 
production of the rj, rj mesons we can see that in the case of the ?/ meson the total cross 
section for the direct production of three pions increases by a factor of about 700 when 
going form rj to rf production threshold. 

The direct three pion production will constitute the major ingredient of the background. 
We have plotted the total cross section of the three pion production using the introduced 
parametrization and compared it to the experimental points with a parametrization of the 
total cross section for the rf meson production. One can see that the cross section for 
the background production is three orders of magnitude larger then the cross section for 
the prodcution of the investigated meson, making the study of the rf — > 7r + 7r~7r° decay 
experimentally challenging. 



E. Parametrization of the proton-proton Final State 
Interaction 



This appendix was prepared on the basis of references [86]. The final state interaction 
between two protons in the exit channel of the reaction pp — > ppX is a well known fact, 
and had achieved a good theoretical description. In order to include this interaction in our 
calculations we had considered three possible models of the pp-FSI enhancement factors 
described in [86]. 

The first model uses the square of the on-shell proton-proton scattering amplitude, 
which for the relative angular momentum I = reads [87]: 

Mpp^p P = — — (E.l) 

o • p 

where p indicates the momentum of both proton in the center of mass system of the colliding 
particles, 5 PP denotes the phase shift, C is the Coulomb penetration factor which includes 
the Coulomb interaction between two protons. This factor can be defined as [88,89]: 

c = ^-r, (E.2) 

where rj c = ^ indicates the relativistic Coulomb parameter, with a being the fine structure 
constant and v the proton velocity in the rest system of other proton. 

The phase shift 5 PP is calculated using the Cini-Fubini-Stanghellini formula, including 
the Wong-Noyes Coulomb corrections [71]: 



C l -p- ctgibpp) + 2-p- Vc - h(r) c ) = -— + ^ ■ b pp ■ p z - ^ P * 2 , (E.3) 

CLpp z i ~r <^pp • p 

where [88]: 



h( Vc ) = -ln( Vc ) - 0.57721 + Vc ■ V - ] , (E.4) 

n=l v L ' 

a pp = —7.8 fm denotes the scattering length and b pp = 2.8 fm is the effective range. The 
parameters P pp = 0.74 fm 3 and Q pp = 3.35 fm 2 are related to the shape of the nuclear 
potential, and were calculated from the pion-nucleon coupling constant [70]. Finally, the 
square of the proton-proton scattering amplitude can be expressed as: 

\Mpp->p P \ 2 = - - - — 2 (E.5) 

As a second possibility, for the pp-FSI parametrization we have approximated the 
enhancement factor by the inverse of the squared Jost function [90]: 

M = M on J- 1 (-p), (E.6) 
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where for a short ranged interaction, by applaying the effective range expansion to order 
p 2 we have: 

Wp 2 + « 2 ) {E7) 



with 



a 



I + a/1 + 2 



(E.8) 



VP 



Equation IE. 71 is true in presence of the long-ranged Coulumb force. 

To include the possible dynamics of the three pion production along with the theoretical 
prediction, we had assumed that the production can proceed via the excitation of A and 
N* resonances according o the reaction chain: 



PP 



AN* — ► ppTTTTTT, 



where A decays into pir and N* into p7rir. A corresponding Feynman diagram could 
then have a form as shown in figure IE. 11 The mechanism had been applied by using the 

7T 
7T 



P- 




P 



7T 



Figure E.l: Feynman diagram of a possible production of three pion by excitation of A and TV* 
resonances. 

relativistic Breite- Wigner form for the probability of the accuring of the m np invariant 
mass: 

2 m 2 -T 2 

(m z — + m z ■ 1 z 

where the V and m denote the width and mass of the A resonance, respectively, and the 
m p7r indicates the invariant mass of the proton-pion subsystem. An analogous formula was 
used for the mr,^ and N* resonance. 



F. Dynamics of the tt + tt tt° production in proton-proton 
interaction 

With a 47r facilities such as WASA, aiming at measurements of decays of rj and 77' pro- 
duced in pp interactions, the understanding of the pp — > pp3ir reaction dynamics becomes 
very important as they constitute a severe background for studies of 77 and rf decays into 
three pions. However, for more than fourty years there were only three experimental data 
points available for the cross section oipp — > ppir + ir~-K and pp — ► pmr + -K + ir~ reactions, all 
coming from bubble chamber experiments [61-63]. Only recently the data base has been 
extended by measurements of the pp — > ppir + ir~'K and pp — > pp7r ir°ir reactions cross 
sections near the threshold by the CELSIUS/ WAS A collaboration [2]. For the remaining 
reactions there are no data in that energy region. As a consequanece the understanding 
of the mechanism of the 3ir production is by far not satisfactory [7]. Therefore, in order 
to gain more information of the dynamics of this process we have compared the empirical 
-jS- energy dependence with results of simulations conducted under various assump- 

tions as regarding the interaction among the final state particles. The direct production 
should proceed by an excitation of one or two baryon resonances followed by the subse- 
quent decays [91]. Therefore we have tested two possibilities assuming for the primary 
production a pure phase space distribution and as a second possibility we considered the 
resonant production via the pp — ► AN* — ► pp3n reaction chain. The result of the simu- 
lations are confronted with the data in figure IF. II Assuming the Watson-Migdal ansatz, 
we factorize the proton-proton FSI and the production amplitude. We have considered 
three possible parametrizations of the pp-FSI enhancement factors which are described in 
detail in reference [86] and in the appendix [E] of this thesis. As a first possibility for the 
proton-proton enhancement factor we use the square of the on-shell proton-proton scat- 
tering amplitude calculated according to the Cini-Fubini-Stanghellini formula including 
the Wong-Noyes Coulomb corrections [70-72]. As a second possibility, we approximated 
the enhancement factor by the inverse of the squared Jost function (Josto) derived by 
Druzhinin et al [90], and finally we used also the inverse of the squared Jost function 
(Jostc), calculated according to the formulae of Goldberger and Watson [92]. Figure [Fj] 
indicates that the excess energy dependence of -^l can be reproduced equally well 

assuming that the phase space is homegeneously populated or assuming that the reaction 
proceeds via resonances: pp — ► AN* — > pp3n. However, it is evident from the figure that 
inclusion of the FSI inhancement factors as derived for the three body final state worsens 
significantly the consistency with the data. 
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Figure F.l: Differential cross section ( ^f| m , ) for the pp — > pp3n production with the invariant 
mass of pions equal to the mass of the rf meson as a function of the excess energy above the rj' 
threshold. Simulations for the pp — ► ppSir reaction are compared to the multimeson production 
extracted from the COSY-11 experiments (see Table |4~T1 and Figure |4~2| . For better visalization 
the simulated values have been arbirarily normalized such that the point at the lowest energy (for 
Q = 1.53 MeV) differ from the data by subsequent powers of ten. The data points are positioned 
in the middle of the figure. The results of simulations were obtained assuming that the process is 
governed by (from top to bottom): 

i) Phase Space © AA © FSIj ostG ~ empty stars 

ii) Phase Space © AA © FSIj ost D - empty triangles 

iii) Phase Space © AA © FSI on -sheii - empty squares 

iv) Phase Space © AA - empty circles 

v) Phase Space - data points 

vi) Phase Space - full circles 

vii) Phase Space © FSI on -sheii - full squares 

viii) Phase Space © FSIj os tD - full triangles 

ix) Phase Space © FSIj os tG - full stars 



G. Description of computer programs used for simulations 

Simulations were performed using the FORTRAN and C language programs with 
additional libraries from the CERNLIB package. Here only a brief description of the 
programs is shown without details. 

All programs were using GENBOD phase space generator subroutines [78], which on 
basis of the total energy in the center of mass of two colliding particles and the masses of 
the outgoing particles, returns the four-momenta of all particles in the exit channel. Each 
generated event is weighted with the WTo parameter which ensures that the particles are 
homogeneneously distributed in the phase space. The output four-momenta are given in 
the center of mass system, thus it is needed to transform them to the laboratory reference 
frame using the Lorentz transformation. 

In most cases the input parameters (beam momentum, total energy in the CM frame, 
number of outgoing particles, mass of particles) for all programs were read from a special 
input file. This method prevented from compiling the code each time, when the values of 
parameters changed. 

In order to check whether the generated events fulfill the detection conditions we have 
implemented constrains of geometrical acceptance of the detector setup. 

The interaction of protons in the final state (FSI) and the production dynamics was 
included according to formulae introduced in appendix [Ej To take into account these 
effects we have weighted each generated event by the sqaure of the proton-proton elastic 
scattering amplitude or a square of the Breit-Wigner amplitude (A 2 ): 

WT = WT -A 2 . (G.l) 

As output of the programs in most cases an . hbook and an ASCI files was produced. The 
first one contains histograms which can be directly ploted by PAW analyzing program [93] 
and the second type of file can be used in further analysis, for example as an imput file to 
another program. 
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"I never did a day's work in my life. It was all fun." 

Thomas Alva Edison (1847 - 1931) 
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